编号 : 2016-0252 


磺 化 聚 共 醚 膜 质子 扩散 的 分 子 动力 学 模拟 研究 
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摘 要 : ”为 了 探究 磺 化 聚 茶 醚 (SPPO) 用 于 质子 交换 膜 燃 料 电池 (PEMFCs) 的 性 质 ， 我 们 研究 了 SPPO 膜 内 部 水 和 水 合 氧 离子 
的 分 子 扩散 特性 。 构 建 含水 量 不 同 的 SPPO 膜 的 分 子 结构 模型 ,采用 分 子 动力 学 (MD) 研 究 了 平衡 态 模型 的 静态 结构 ， 发 现 元 胞 
的 密度 随 着 含水 量 的 增加 而 降低 ， 同 时 通过 分 析 径 向 分 布 函数 (RDF)， 发 现在 SPPO 链 上 的 磷酸 根 上 的 硫 原子 周围 有 水 分 子 聚 
集 ， 起 到 吸水 和 传导 的 作用 ， 水 合 氢 离 子 也 更 趋向 于 聚集 在 磺 酸 根 附近 。 通 过 对 分 子 运动 轨迹 的 分 析 ， 计 算 了 水 和 水 合 氨 离 子 
在 SPPO 膜 内 的 扩散 系数 D 和 电导 率 c， 研 究 结果 表 明 ，SPPO 膜 的 电导 率 随 着 膜 内 水 含量 的 增加 而 增加 ， 得 到 了 与 实验 符合 
较 好 的 结果 。 
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中 图 分 类 号 : TK124 文件 标识 码 :: A 


Molecular Dynamics Simulation Study of Proton Diffusion in Polymer 


Electrolyte Membranes Based on Sulfonated Poly (phenylene oxide? 


Han Xuhui, Zhang Peng, Ma Yue, Li Jia, Wang Yuxiang, Zhou Qiong, Zhou Guanggang, Lu Guiwu 
(China university of petroleum (Beijing), college of science; Beijing 102249) 

Abstract To explore the possibility of sulfonated poly (phenylene oxide) (SPPO) membrane for proton exchange 
membrane fuel cells(PEMFC), molecular dynamics(MD) simulation technique was employed to investigate the 
molecular diffusion characteristics of water and hydrogen ion in SPPO membrane of different water content. 
Through analysis the structure of the fully equilibrated amorphous polymeric cell, we find that the density decreased 
with the increase of water contents. According to the radial distribution function (RDF) of cell, it was found that the 
water and hydrogen ion molecules tend to gathered in the sulfonate group of SPPO chains which made the sulfonate 
group take the conductive roles. Through the dynamic property analysis, we calculate the diffusion coefficient of 
water and hydrogen ion in the SPPO membrane, and the conductivity of SPPO membrane. With the increase of water 
content, the diffusion coefficient of water and hydrogen ion and the conductivity also increases. The results coincide 
well with the experiments. 
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TESTES 高 成 本 ， 是 阻碍 其 应 用 于 直接 甲醇 燃料 电池 


燃料 电池 具有 安全 、 高 效 、 无 污染 、 适 用 广 、 
无 噪声 等 特点 ， 已 成 为 当今 世界 能 源 领域 的 研究 热 
点 ! 了 。 在 几 类 重要 的 燃料 电池 分 类 中 ， 聚 合 物 燃 料 
电池 (PEMFCs) 受到 了 广泛 的 关注 口 。 

质子 交换 膜 是 一 种 由 高 分 子 材料 制 成 的 含有 离 
子 交 换 基 团 的 一 种 阳离子 交换 薄膜 ， 在 PEMFCs 中 
起 着 质子 传导 、 隔 离 燃 料 如 氧气 和 甲醇 与 氧气 和 绝 
缘 电 子 的 作用 , 它 是 PEMFCs 研究 的 重点 和 关键 1。 
目前 ,只 有 全 气 磺 酸 型 质子 交换 膜 (商品 名 为 Nafion) 
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(DMFCs) 的 主要 原因 1。 磺 化 聚 茶 醚 膜 (SPPO) 因 
具有 和 较 独 特 的 性 质 而 引起 了 广泛 的 关注 总。 以 往 对 
SPPO 膜 的 研究 主要 集中 在 作为 气体 分 离 膜 用 于 反 
ARR, HA, BARDA SPPO 膜 作 为 电解 
质 材料 用 于 燃料 电池 的 研究 报道 号 -六 。 

近 几 年 ， 对 于 质子 交换 膜 的 分 子 动力 学 模拟 研 
究 已 有 大 量 的 报道 后 。 从 理论 上 模拟 及 预测 SPPO 
膜 性 质 ， 显 得 很 有 必要 504。 本 文 利用 Materials 
studio(MS) 软 件 平台 构建 了 SPPO 膜 不 同 含水 量 的 结 
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构 模 型 ， 模 拟 了 水 分 子 和 水 合 氧 离子 在 膜 内 的 扩散 
行为 ,计算 了 对 应 的 电导 率 。 
1 计算 模型 与 模拟 方法 

为 了 研究 水 和 水 合 氧 离子 在 SPPO 膜 内 的 传导 
特性 ， 构 建 了 磺 化 度 为 50%， 聚 合 度 为 10 的 SPPO 
链 ， 如 图 1 所 示 。 其 中 图 1(a) 和 (b) 分 别 为 未 磺 化 和 
磺 化 的 SPPO 单 体 结构 模型 ， 图 1(c) 为 SPPO 链 的 构 
型 。 


图 1 未 磺 化 (a) 和 磺 化 (b) 的 SPPO 单 体 的 模型 和 SPPO 链 (c) 
的 分 子 模型 

Fig. 1 Chemical structures of non-sulfonated (a) and sulfonated 

(b) SPPO monomers used in the simulation. and SPPO chain (c) 


molecular model 


利用 Amorphous cell 模块 构建 SPPO 膜 的 质子 传 
导体 系 ， 图 2 为 水 和 磺 酸 比 A-5. 8. 11. 14 (4= 水 
分 子 个 数 / 磺 酸根 个 数 ) 的 SPPO 膜 的 分 子 模型 ， 具 
体 元 胞 的 参数 如 表 1 所 示 。 
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图 2 4=5(a)、8(b)、11(c)、14(d) 的 元 胞 结构 
Fig2 The cell structure for À-5(a). 8(b). 11(c)、14(d) 


表 1 元 胞 参数 


Table 1 parameters of cell 
名 称 À 水 分 子 个 数 “” CKAO€XBTTA 
SPPO-5 5 200 50 
SPPO-8 8 350 50 
SPPO-11 11 500 50 
SPPO-14 14 650 50 


利用 MS 软件 中 的 Forcite 模块 进行 模型 优化 和 
分 子 动力 学 模拟 , 所 有 动力 学 模拟 均 选 用 COMPASS 
力 场 。COMPASS 是 以 从 头 计 算 为 依据 的 力 场 ， 采 用 
从 头 计算 和 经 验方 法 相 结合 的 方式 来 计算 力 场 中 的 
参数 。 它 是 第 一 个 把 有 机 分 子 体系 的 力 场 与 无 机 分 
子 体系 的 力 场 统一 的 分 子 力 场 ， 能 够 模拟 有 机 和 无 
机 小 分 子 、 高 分 子 、 一 些 金属 离子 、 金 属 氧化 物 与 
金属 71， 

初始 模型 构建 完成 后 ， 由 于 体系 中 的 分 子 初始 
构象 并 不 是 能 量 最 低 的 构象 ， 因 此 需要 利用 退火 过 
程 对 模型 进行 优化 处 理 。 采 用 最 速 下 降 法 对 元 胞 结 
构 的 能 量 和 几何 结构 进行 优化 ， 通 过 退火 模块 模拟 
退火 过 程 找到 模型 能 量 最 低 构 象 。 最 后 ， 在 正则 系 
综 下 ， 设 置 体系 温度 为 298K, 动力 学 时 间 1.5ns, 2 
fT SPPO 膜 质 子 传导 的 分 子 动力 学 模拟 。 本 文采 用 的 
时 间 模 拟 步 长 均 为 1fs。 

通过 以 上 步骤 得 到 不 同 含水 量 的 模型 的 动力 学 
运动 轨迹 ， 对 其 最 终 的 模型 结构 和 运动 轨迹 进行 分 
析 ， 得 到 元 胞 密度 、 体 积 和 均 方位 移 。 


2 结果 与 讨论 


2.1 元 胞 密度 分 析 

对 于 非 晶 态 元 胞 的 动力 学 模拟 来 说 ， 最 重要 的 
是 要 保证 模拟 元 胞 要 处 于 完全 平衡 状态 ， 以 确保 水 
分 子 和 水 合 氢 离子 能 够 充分 地 扩散 到 处 于 平衡 态 
SPPO 膜 中 。 为 此 引入 温度 、 总 能 量 和 密度 等 参数 ， 
来 观察 体系 是 否 已 经 处 于 平衡 态 。 图 3 为 体系 的 温 
度 T 和 总 能 量 E 随时 间 1 变化 的 波动 曲线 ,图 中 曲线 
波动 很 小 即 体 系 温度 和 总 能 量 基 本 保持 不 变 ， 表 2 
为 模拟 的 元 胞 的 体积 和 密度 ， 图 3 和 表 2 反映 出 体 
系 已 经 处 于 平衡 状态 。 
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3 模拟 体系 温度 和 总 能 量 随时 间 变 化 曲线 
Fig.3 Temperature and total energy of the dynamic 


simulation system 


42 模拟 元 胞 的 密度 和 体积 


Table 2 density and volume of cell 


膜 名 称 密度 /g/cm 体积 /nm 
SPPO -5 1.24 27.50 
SPPO -8 1.21 32.00 
SPPO-11 1.17 36.90 
SPPO-14 1.14 40.76 


2.2 静态 结构 分 析 
图 4(a)、(b)、(c) 为 三 5 的 SPPO-5 膜 的 径 向 分 布 

函数 (Radial distribution function，RDF)。 将 人 磺 酸 
根 上 的 硫 原 子 标记 为 S， 水 分 子 上 的 氧 原子 标记 为 
Ow， 水 合 氨 离 子 上 的 氧 原 子 标 记 为 O， 甲 基 上 的 碳 
原子 标记 为 Cm KA ERRETEN Coo ER 
4(a) 中 的 两 了 个 曲线 我 们 发 现 S-Oy 曲线 和 S-O,, 曲线 
有 明显 的 峰值 ，S-Ou 曲线 比 S-O 曲线 的 峰值 高 。 这 
是 因为 水 合 氢 离子 所 带 的 正 电 与 带 负 电 磺 酸根 基 团 
相互 吸引 导致 水 合 氢 离 子 更 多 的 在 人 磺 酸 根 附近 存 
在 。 将 图 4(a) 中 的 曲线 与 图 4(b) 比 较 我 们 发 现 ， 在 
SPPO 链 上 水 合 氢 离子 更 趋同 于 在 磺 酸 根 周 围 ， 而 不 
是 甲烷 和 葵 环 等 疏水 基 团 。 图 4(c) 中 我 们 看 到 S-S、 
Orn-On 的 截断 半径 在 0.4 nm 以 后 ， 而 图 4(a) 中 S-O,. 
S-O, 的 截断 半径 在 0.4 am， 这 也 能 说 明 磺 酸 根 和 水 
合 氢 离 子 相 互 吸 引 ， 相 对 于 自身 水 合 氢 离 子 更 趋 癌 
聚集 在 磺 酸 根 附近 。 对 比 4=8、11、14 的 SPPO 膜 
的 RDF 曲线 也 能 发 现 相 似 的 规律 。 
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图 4 (a) SPPO-5 膜 中 S-O,、S-O。 径 向 分 布 函数 ;(b) SPPO-5 
膜 中 Cu-O、 Cn-On 径 向 分 布 函数 ; (c) SPPO-5 膜 中 S-S、Oh-Oh 
径 向 分 布 函数 
Fig.4 (a) RDF of S-O,. S-O,, in SPPO membrane;(b) RDF of 
Cp-On~ Cm-On in SPPO membrane ; (c) RDF of S-S, O,-O, in 

SPPO membrane 


2.3 水 和 水 合 所 离子 在 膜 中 的 扩散 系数 计算 
初始 模拟 体系 经 过 1.5ns 的 动力 学 模拟 后 ， 取 动 
力学 过 程 100ps-1200ps 时 间 段 ， 利 用 分 析 工 具 得 到 
如 图 5 所 示 的 均 方 位 移 图 。 由 于 模型 中 水 与 水 合 氢 
原子 的 体积 均 很 小 ， 所 以 采用 记录 水 和 水 合 氧 原子 


中 氧 原 子 的 位 置 数据 来 计算 扩散 系数 。 将 水 中 的 氧 
(O,2 和 水 合 氨 离 子 中 的 氧 COD 分 别 标识 。 
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5 (a) 不 同 舍 水量 水 分 子 的 MSD; (b) 不 同 含水 量 水 合 氧 
离子 的 MSD 


Fig. 5 (a) MSD of water at different water content; (b) MSD of 


hydronium ions at different water content 


扩散 系数 DD 

a 
D=— 1 
(1) 


式 CD 中 : a 为 对 数据 进行 拟 合 得 到 均 方位 移 (MSD) 
与 时 间 的 一 次 函数 的 曲线 的 斜率 ， 将 a 代入 式 (1) 中 
得 到 扩散 系数 D。 

ILI 


2 
B Dne 


° =T Q) 


A OD 中 ;为 玻 尔 兹 曼 常 数 ，n 为 元 胞 中 水 合 氧 离 
子 的 粒子 个 数 ，7 为 绝对 温度 ，e 为 元 电荷 电量 。 将 
DD 代入 式 (2)， 可 计算 得 到 体系 的 电导 率 。 

表 3 为 最 终 得 到 的 水 分 子 和 水 合 氢 离子 的 扩散 
系数 及 电导 率 ， 从 表 3 中 可 以 看 出 随 着 含水 量 的 增 
加 ， 水 和 水 合 氢 离子 的 扩散 系数 随 之 增加 。 反 映 了 
水 含量 越 多 ， 膜 就 越 湿 洞 ， 质 子 的 传导 率 就 越 大 的 


客观 事实 。 将 水 合 氢 离子 的 扩散 系数 代入 到 式 (2)， 
计算 得 到 相应 的 电导 率 如 表 3 所 示 。 由 表 3 容易 看 
H, 随 着 扩散 系数 的 增 大 , SPPO 膜 的 电导 率 就 越 大 ， 
变化 范围 从 0.0056-0.0382 S/cm. 至 春 丽 等 (1 对 SPPO 
膜 进行 的 电导 率 测试 , 测试 磺 化 度 40.1% 的 SPPO 膜 
时 ， 其 电导 率 为 0.0116 Secm， 与 我 们 的 模拟 结果 基 
本 符合 。 


表 3 水 分 子 、 水 合 氨 离 子 的 扩散 系数 和 电导 率 
Table 3 Diffusion coefficient of water and 


hydronium ions and conductivity 


膜 名 称 4 D, D, o/ S/cm 
SPPO -5 5 0.0133 0.0050 0.0056 
SPPO -8 8 0.0433 0.0150 0.0127 
SPPO-11 11 0.0783 0.0250 0.0185 
SPPO-14 14 0.1067 0.0500 0.0382 


3 Ae 


对 SPPO 膜 的 元 胞 进行 了 分 子 动力 学 模拟 , 研究 
结果 表明 ，SPPO 膜 的 密度 随 水 含量 的 增加 而 下 降 ， 
磺 酸 根 承 担 质子 传导 的 作用 ， 水 合 氢 离子 相对 更 容 

聚集 到 磺 酸 根 附 近 。 研 究 发 现 水 分 子 和 水 合 氧 离 
子 的 扩散 系数 随 着 水 含量 的 增 大 而 变 大 ， 其 电导 率 
也 是 随 着 水 含量 的 增加 而 增加 ， 其 规律 与 实验 基本 
相符 。 
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